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James  F.  Driscoll 

Assistant  Professor,  Department  of  Aerospace  Engineering 
The  University  of  Michigan,  Ann  Arbor,  Michigan  48109 


A.  . ABSTRACT 

< The  one  year  research  program  described  here  was  a direct  extension  of 
work  completed  during  the  1977  AFOSR/ASEE  Summer  Faculty  Research  program 
at  AFRPL.  During  that  program  the  author  and  Dr.  David  Mann,  AFRPL,  performed 
the  first  in  situ  size  measurements  of  submicron  alumina  (A1|0^)  particles 
that  have  been  published  to  date.  The  new  laser  diagnostic  technique  used. 
Diffusion  Broadening  Spectroscopy,  represents  a new  advance  in  the  state- 
of-the-art  of  submicron  particle  sizing. 

The  goal  of  the  research  described  herein  was  to  develop  a data  reduction 
procedure  to  extend  the  applicability  of  the  technique,  so  that  size  distribu- 
tion half  width  can  be  measured,  in  addition  to  particle  mean  diameter. 

The  computer  routine  was  developed  to  determine  particle  mean  diameter 
and  size  distribution  half  width,  given  light  scattering  data  obtained  at  two 
scattering  angles.  Sensitivity  of  the  results  was  investigated  and  feasibility 


3.  Multiple  angle  scattering,  using  three  or  more  detectors,  as  a 
method  of  measuring  particle  size  distribution  was  considered. 
Analysis  was  performed  to  indicate  the  feasibility  of  obtaining 
particle  size  distribution  directly,  without  assuming  any  func- 
tional form  of  the  distribution. 


C.  STATUS  OF  THE  RESEARCH  EFFORT 

Progress  toward  achieving  the  three  major  research  objectives  is  now 
outlined. 

1.  Aerodynamic  drag  on  submicron  particulates  was  investigated  for 
gas  pressures  between  10"^  to  1 atmosphere.  Current  theories 
describing  free  molecular  flow  drag  forces,  transition  flow,  and 
continuum  flow  forces  were  assessed.  Of  the  numerous  theories, 
the  most  realistic  is  one  proposed  by  N.A.  Fuchs  which  agrees 
with  experiments  performed  by  Millikan  to  within  2%.  Particle 
drag  in  the  free  molecular,  transition,  and  continuum  ranges  can 
be  expressed  as: 

Drag  = 6 irydU  [(1  + .42  t/d)-1  + 1.67  i/d]"1 
Some  values  of  particle  drag  are  given  in  Fig.  1. 

2.  A distribution  of  particle  sizes  was  considered  and  the  net 
diffusion  was  investigated.  A first  step  was  to  assume  a func- 
tional form  of  the  distribution  (i.e.  a log-normal  distribution 
shape)  and  to  determine  if  the  geometric  mean  diameter  and 
geometric  standard  deviation  could  be  deduced  using  the  DBS 
measuring  technique.  The  results  show  such  measurements  are 
possible,  if  the  scattering  angle  is  optimized  and  the  laser 
wavelength  is  chosen  to  be  at  least  ten  times  the  particle 
diameter  (Mie  scattering  regime).  Details  of  how  to  measure  the 
geometric  mean  diameter  and  standard  deviation  of  the  aerosol  are 
described  in  the  attached  publication.  The  major  portion  of  the 
work  done  on  this  program  is  outlined  in  this  publication.  It  was 
necessary  to  develop  an  extensive  computer  program  to  account  for  a 
polydisperse  aerosol  and  low  pressure  effects  as  mentioned  above  in 
Item  C.l . 
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The  kinetics  of  transport  processes  for  intermediate  Kn 


Figure  1.  Particle  Drag  at  Lnw  Pressures 

(Fuks,  Highly  Disbursed  Aerosols) 


3.  Multiple  angle  scattering  was  considered  in  order  to  determine  the 
particle  size  distribution  curve  without  assuming  a log-normal 
shape  to  the  curve.  Instead  of  a two  parameter  log-normal  curve,  a 
more  general  ten  parameter  curve  was  considered,  requiring  inde- 
pendent measurements  at  ten  different  scattering  angles.  Analysis 
showed  that  the  uncertainty  of  results  obtained  using  this  method 
would  be  too  high  for  practical  application  unless  extreme  care  was 
taken  to  optimize  the  accuracy  of  each  of  the  ten  measurements. 

This  did  not  appear  feasible  within  the  limitations  of  most  university 
research  programs. 


0.  PUBLICATIONS 

"Submicron  Particle  Size  Measurements  in  an  Acetyl ene/Oxygen  Flame," 

J.F.  Driscoll,  D.M.  Mann,  Third  Int'l  Workshop  on  Laser  Velocimetry, 
Purdue,  1978. 

"Submicron  Particle  Size  Measurements  in  an  Acetylene-Oxygen  Flame," 

J.F.  Driscoll,  D.M.  Mann,  W.K.  McGregor,  accepted  for  publication, 
Combustion  Science  and  Technology. 

E.  PERSONNEL 

James  F.  Driscoll,  Assistant  Professor 

F.  OTHER  PROGRESS: 

Equipment  has  been  purchased  using  University  funds  to  allow  this 
research  effort  to  be  continued.  A submicron  particle  generator  system 
and  an  argon  ion  laser  have  been  purchased.  A proposal  to  AFOSR  for  continu- 
ing research  in  this  area  is  planned. 
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ABSTRACT 

Diffusion  broadening  spectroscopy  has  been  used  to  make 
in  situ  size  measurements  of  submicron  aerosol  particles  in 
an  acetylene-oxygen  flame.  Making  use  of  the  spectjral  broad- 
ening of  scattered  laser  light  due  to  random  particle  motion, 
this  optical  technique  has  been  demonstrated  to  be  suitable 
for  application  in  unsteady  high  temperature  environments  of 
2,200°K  in  flowing  gases.  Particle  sizes  were  found  to  vary 
from  40  to  250  nm  as  the  height  above  the  burner  varied  from 
0.5  to  3.0  cm.  Little  variation  was  found  as  a result  of 
equivalence  ratio  changes  from  2.5  to  5.0.  Sizes  measured 
from  electron  microscope  photographs  of  soot  particles  col- 
lected from  the  flame  show  reasonable  agreement  with  those 
determined  from  scattering  measurements.  A method  to  mea- 
sure geometric  mean  diameter  and  standard  deviation  of  the 
particle  size  distribution  in  a polydisperse  aerosol  is  also 
presented,  and  experimental  results  described.  Modifications 
necessary  to  account  for  polydispersity  and  noncontinuum  flow 
effects  are  discussed. 
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INTRODUCTION 

The  study  of  particle  formation  in  flames  and  other  conden- 
sation processes  has  been  limited  by  the  lack  of  adequate  meth- 
ods to  measure  particle  size.  Several  optical  methods  have  been 
proposed  to  determine  particle  size  in  situ,  and  thereby  elimi- 
nate the  problems  inherent  to  conventional  sampling  probes,  which 
can  perturb  the  particle  agglomeration  and  reaction  process. 
However,  flame  particulates  present  special  problems  since  they 
generally  occur  in  the  submicron  size  range  and  are  usually  of 
unknown  composition.  For  this  reason,  the  use  of  sizing  diagnos- 
tics such  as  holography  and  fringe  visibility  interferometry 
(which  are  limited  to  approximately  5 ym)  and  Mie  scattering 
(which  depends  on  particle  composition)  has  been  limited  in  flames. 
Until  recently,  no  single  technique  has  afforded  the  capability  of 
in  situ  sizing  of  submicron  particulates  of  unknown  composition. 

A relatively  new  optical  technique  that  can  provide  this 
capability  was  demonstrated  by  Cummins  et  al.  (1964)  and 
Foord  et  al.  (1970)  in  liquids,  and  later  by  Hinds  and  Reist 
(1972)  in  a gas.  This  technique,  which  will  be  referred  to  as 
Diffusion  Broadening  Spectroscopy  (DBS) , relies  on  physical 
principles  which  were  outlined  by  Cummins  et  al.  (1964) , Benedek 
(1969)  and  Edwards  et  al.  (1971) . Particle  size  is  determined 
directly  by  measuring  the  spectral  broadening  of  scattered  laser 
light.  Broadening  of  the  observed  optical  spectrum  is  a result 
of  two  factors:  the  vigorous  diffusion  of  particles  due  to  their 
Brownian  motion  and  the  Doppler  effect.  Therefore,  the  particle 
diffusion  coefficient  can  be  determined  directly  from  a measure- 
ment of  the  halfwidth  of  the  optical  spectrum.  Particle  size  is 
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then  inferred  using  a modified  Stokes-Einstein  relation  that 
relates  the  diffusion  coefficient  and  the  particle  diameter. 

Several  calibration  studies  have  verified  that  accurate 
results  can  be  obtained  using  the  DBS  technique. 

Cummins  et  al.  (1964)  used  dilute  solutions  of  monodisperse 
.088  iim  polystyrene  molecules  and  Hinds  and  Reist  (1972)  used 
airborne  aerosols  of  monodisperse  polystyrene  particles  in 
the  .557  to  1.10  micron  range.  Penner  et  al . (1976)  modified 
the  technique  for  use  in  flowing  gases  and  presented  results 
obtained  in  an  ethylene-oxygen  flame. 

In  the  work  described  herein,  particle  size  measurements 
have  been  performed  at  various  locations  in  an  acetylene- 
oxygen  flame,  using  Diffusion  Broadening  Spectroscopy.  Results 
are  compared  with  particle  sizes  determined  from  electron 
micrographs  of  particles  collected  from  the  flame.  Methods 
to  optimize  signal  to  noise  ratio  in  combustion  processes  and 
thereby  enhance  the  range  of  size  measurements  possible  have 
been  identified.  In  order  to  size  flame  particulates,  two 
modifications  to  the  DBS  technique  were  required.  For  our 
flame  conditions,  it  was  necessary  to  account  for  aerosol 
polydispersity  and  free  molecular  flow  effects. 


THEORETICAL  CONSIDERATIONS 


The  fundamental  physics  underlying  Diffusion  Broadening 
Spectroscopy  is  described  below.  This  treatment  builds  upon 
theory  presented  in  Hinds  and  Reist  (1972).  A particle  j, 
located  at  position  vector  r in  the  illuminating  field  of 
laser  radiation  cf  frequency  <i)o,  scatters  light  to  a detector 
located  at  a scattering  angle  0 from  the  incident  beam.  The 
scattered  electric  field  at  the  detector  is: 

E^  (t)  = Eq  exp  (i(k-r(t)  - mot) ] (1) 

where  Eq  is  the  magnitude  of  the  field  at  the 
detector,  and  the  scattering  wave  vector  k is: 

k = [(2tt/Xo)  (e±  - es)  1 (2) 

Vectors  e.  and  e are  unit  vectors  in  the  incident  and  scattered 
i s 

wave  directions,  respectively.  The  magnitude  of  k is 

K = [ (2tt/A0)  (2  sin  0/2)]  (3) 

Particle  motion  results  in  a frequency  shift  of  the 
scattered  light  which  is  represented  by  the  time  variation  of 
the  phase  factor  k*r  in  Eq.  (1) . For  the  case  of  uniform 
particle  velocity  v,  the  phase  factor  becomes 

k*r  = k-vt  + k-7o  (4) 

resulting  in  a constant  frequency  shift  u>D  given  by 
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(5) 


u>D  =|k*v|=  (2Tr/AQ)  v (eg  - e.^) 


which  is  the  familiar  Doppler  shift  relation  used  in  laser 
Doppler  velocimetry.  The  vector  rQ  simply  defines  the 
coordinate  system  location. 

To  assess  the  effects  of  particle  motion  on  the  power 
spectrum  of  scattered  light,  the  field  correlation  function 
is  first  computed: 


E*(t)  E (t+x) 
s s 


= E E E £ E 3 
l j ° o 


exp  [ik  * (r j (t)  - r..  (t  + t))] 


where  E is  the  sum  of  individual  contributions:  . 
s 


E (t)  = E E ^ (t) 
S j=l  S 


Under  the  assumption  that  the  scattering  from  each  particle  is 
independent,  that  the  illumination  is  uniform  and  that  the 
position  of  each  particle  depends  only  on  the  time  interval  t: 


^(t)  = E*(t)  Eg  (t  + t)  - N exp  [i  k • Ar  ] 


Ar j = r j (t)  - r j (t  + T ) 


The  average  over  Ar  is  given  by: 


exp  (i  k *Ar) 


■ I 


f(Ar,t)  exp  (i  k *Ar)  dAr 


■'AmI  'X 'Pf  . 


where  f(Ar,t)  is  the  distribution  of  displacements  Ar  during 
the  time  interval  t • Furthermore,  the  intensity  correlation 
function,  which  is  of  physical  interest,  is  related  to  the 
field  correlation  function  R (t)  by: 

1j 


Rj.(x) 


(11) 


At  this  point  it  is  necessary  to  incorporate  the  physics 
of  particle  movement.  For  a particle  undergoing  Brownian 
motion,  the  probability  of  a displacement  Ar  in  a time  t is: 


f(Ar,x)  = (4ttDt)  'L//2  exp  (-  [ Ar  1 2/4  Dx) 


(12) 


where  D is  the  particle  diffusion  coefficient.  The  resulting 
optical  spectrum  can  now  be  obtained  from  the  intensity  corre- 
lation function  by  applying  the  Wiener-Khinchine  theorem: 


oo 


Pi(u>)  = 


(1/n) 


Ra  (x) 


exp  (iaix)  dx 


(13) 
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Combining  Eqs.  (8)  through  (13)  gives  the  form  of  the  optical 
spectrum: 


P^w) 


2 2 

E K D/tt 
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(K2D)2  + („,  - u>  ) 2 


(14) 


which  is  a Lorentzian  centered  at  the  laser 


frequency  wq. 


The  halfwidth  of  optical  spectrum  given  by  Eq.  (14)  is 
measured  using  homodyne  detection.  For  the  case  of  many 
particles  radiating  simultaneously  in  the  focal  volume,  all 
frequency  components  in  the  broadened  optical  spectrum  are 
simultaneously  incident  on  the  photodetector.  The  beating  of 
each  discrete  frequency  interval  of  the  optical  spectrum  with 
every  other  discrete  frequency  interval  results  in  a spectral 
profile  centered  at  zero  frequency,  called  the  homodyne  spec- 
trum, which  is  a Lorentzian  function  for  a monodisperse, 
stationary  aerosol: 


Sj  (w) 


I2  2K2D 

11  (2K2D)2  + co2 


(15) 


The  homodyne  spectrum  described  by  Eq.  (15)  has  a halfwidth 
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at  half  height  of  2K  D radians/sec,  corresponding  to  a half- 
width of  K2D/tt  Hertz. 

The  particle  diameter  d in  a monodisperse  aerosol  is 
related  to  the  diffusion  coefficient  D by  the  modified 
Stokes-Einstein  relationship: 

D=3^C|ll/d)  (16> 

where  the  temperature,  viscosity  and  mean  free  path  of  the  gas 
are  designated  T,  p and  i,  respectively.  The  effect  of  mean 
free  path  £.  in  Eq.  (16)  becomes  important  in  the  transition 
flow  regime  0.1  < Jt/d  < 10  in  which  neither  continuum  nor  free 
molecular  drag  models  apply;  the  slip  correction  factor  C(Jt/d) 
is  given  by  Fuchs  (1970). 


Cowhi  ninq  Eqs . (3),  H 5)  , a.', cl  (16)  , yields  a relation 


for  the  particle  diameter  in  a monodisperse  aerosol: 


d 
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kT  . 2/  0\  C 

2 Sin  ^ 2 HW 
X \i 


(17) 


where  the  halfwidth  HW  at  half  height  is  measured  in  Hertz. 

To  apply  Eq.  (17) , operating  parameters  must  be  selected  to 
minimize  the  effects  of  other  spectral  broadening  mechanisms. 
One  of  the  most  serious  considerations  is  the  minimization  of 
time  of  flight  broadening  (Wang,  1970) , which  occurs  due  to 
the  finite  duration  of  the  scattered  signal  from  each  particle 
as  it  traverses  the  laser  beam.  The  power  spectrum  of  a sine 
wave  of  finite  duration  t is  not  a delta  function,  but  has  a 
spectral  halfwidth  inversely  proportional  to  x , which  corres- 
ponds to  the  particle  transit  time.  The  Gaussian  spectrum 
due  to  time  of  flight  broadening,  therefore,  has  a halfwidth 
equal  to  U/2a,  where  U is  the  flow  velocity,  and  a is  the 
streamwise  Gaussian  halfwidth  of  the  optical  sampling  volume. 
To  apply  the  DBS  technique  to  high  velocity  flow  fields, 
the  local  laser  beam  diameter  a must  be  increased  (thereby 
sacrificing  spatial  resolution)  to  a value  such  that: 


a >>  U/  (4K2D) 


(18) 


Turbulence  and  mean  flow  gradients,  can  also  broaden  the 
homodyne  spectrum,  but  their  effects  can  be  minimized  by 
insuring  that  the  Doppler  shift  due  to  the  maximum  difference 

in  particle  velocities  that  occur  simultaneously  in  the  detec- 

. 2 
tion  volume  is  much  less  than  2l<  D.  Such  velocity  differences 

can  be  reduced  by  reducing  the  beam  diameter,  if  necessary. 


A basic  requirement  of  the  laser  light  source  is  that  the 
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laser  line  width  be  significantly  less  than  the  halfwidth  of 
the  homodyne  spectrum.  This  requirement  is  met  by  conventional 
helium-neon  and  argon  lasers,  for  which  the  linewidth  of  each 
longitudinal  mode  in  the  laser  output  is  less  than  10  hertz. 

Due  to  thermal  instabilities,  the  centerline  frequency  of  each 
longitudinal  mode  fluctuates  and  results  in  a time  averaged 
linewidth  of  each  mode  of  approximately  10  KHz  for  helium 
neon  and  10  MHz  for  argon-ion  lasers.  These  fluctuations  have 
no  effect  on  the  homodyne  intensity  spectrum,  which  is  a 
result  of  the  beating  of  different  frequency  components  which 
must  occur  simultaneously.  Furthermore,  multimode  laser 
operation  has  no  affect  on  the  halfwidth  of  the  homodyne 
spectrum  if  this  halfwidth  is  less  than  the  laser  mode  spacing, 
which  is  typically  200  MHz. 

EXPERIMENTAL  CONDITIONS 

One  inherent  advantage  of  the  use  of  the  DBS  technique 
is  the  simplicity  of  the  experimental  system,  making  it 
especially  attractive  for  use  in  hostile  environments. 

Scattered  radiation  from  a single  laser  beam  is  collected  at 
a selected  scattering  angle  and  measured  using  a photomulti- 
plier tube.  The  power  spectrum  of  phototube  current  fluctua- 
tion is  then  measured  using  a conventional  spectrum  analyzer 
or  by  the  photon  correlation  method  (Cummins  and  Pike,  1973)  . 


A schematic  of  the  experimental  configuration  used  in 
this  study  is  shown  in  Fig.  1.  A 28  milliwatt  helium  neon 
laser  was  used  with  a beamsplitter  to  provide  two  beams  of 
equal  intensity  which  were  focussed  to  intersect  at  the 
sampling  location.  This  arrangement  allowed  scattered 
radiation  to  be  observed  at  either  a 25°  or  60°  angle  from 
the  incident  beam  using  fixed  detection  optics,  by  alternately 
blocking  one  of  the  beams.  Local  beam  diameter  was  1 mm.  A 
collecting  lens  of  15.2  cm  focal  length  was  used  to  focus  the 
enlarged  (3X)  image  of  the  focal  volume  on  the  surface  of  an 
RCA  4832  photomultiplier  tube.  An  interference  filter  cen- 
tered at  632.8  nm  with  a measured  halfwidth  of  1.6  nm  was 
used  to  filter  out  flame  emission.  Photomultiplier  current 
was  passed  to  ground  through  a .1  MH  resistor,  providing  a 
phototube  frequency  response  of  50  KHz. 

Scattered  laser  radiation  from  flame  particulates  was 
visually  observed  through  a pinhole  drilled  in  the  back  wall 
of  the  phototube  mount,  and  proper  focus  was  obtained  when 
the  entire  collection  lens  (i.e.  detection  solid  angle)  was 
filled  with  scattered  laser  light.  The  ratio  of  signal  to 
flame  emission  was  measured  to  exceed  25  and  the  measured 
ratio  of  signal  to  shot  noise  was  greater  than  40. 

A flat  flame  laboratory  burner  was  used  in  this  study 
which  consisted  of  a core  of  bundled  copper  rods  between 
which  flowed  premixed  acetylene,  oxygen  and  a nitrogen  dilu- 
tant introduced  to  vary  flame  temperature  and  thus  flame  speed 


The  3.18  cm  diameter  burner  surface  was  recessed  and  surrounded 


by  an  annular  flow  of  nitrogen,  which  considerably  aided  flame 
stability.  The  entire  apparatus  was  water  cooled.  Unburned  gas 
velocity  was  determined  using  a flowmeter  to  be  20  cm/sec  at  a 
fuel  air  equivalence  ratio  <f>  = 2.5. 

The  measured  flame  emission  was  relatively  steady,  with 
+ 10%  variations  at  frequencies  of  5-15  Hertz.  The  intensity 
of  scattered  laser  light  was  unsteady  due  to  time  varying 
aerosol  concentrations.  However,  unsteady  aerosol  conditions 
pose  no  fundamental  limitations  to  the  sizing  technique  since 
all  signal  fluctuations  of  frequency  less  than  1 KHz  were 
later  removed  by  filtering  during  the  data  reduction. 

Fluctuations  in  the  phototube  current  were  recorded  on 
analog  tape  using  an  AMPEX  1200  FM  recording  system  and 
analyzed  using  a digital  spectrum  analyzer  (Spectral  Dynamics 
360-35) . This  unique  instrument  displays  the  time  averaged 
power  spectrum  as  recorded  data  is  continuously  digitized  and 
Fourier  transformed.  Data  was  first  high  pass  filtered  at 
1 Khz  and  digitized  at  sampling  rates  of  16  - 32  KHz.  Each 
computed  spectrum  was  determined  from  2.5  x 10^  data  points. 

PARTICLE  SIZE  DETERMINATION 

The  size  of  particles  in  the  acetylene-oxygen  flame  was 
determined  from  measured  homodyne  spectra  using  Eq.  (17).  An 
example  of  a filtered  spectrum  is  shown  in  Fig.  2.  The  half- 
width (HW)  of  the  Lorentzian  homodyne  spectrum  is  related  to 
the  filtered  spectrum  by  the  relation: 
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Hw2  ■ fHM2  - 2£o2  <18» 

where  f is  the  frequency  in  Hertz  at  which  the  intensity 
HM 

spectrum  has  decreased  to  half  the  maximum  measured  value, 
and  f is  the  high  pass  filter  frequency,  which  in  this 
case  was  1 KHz. 

Flame  temperature  and  concentration  of  26  combustion 

species  were  calculated  using  the  AFRPL-ISP  computer  routine, 

which  utilizes  updated  JANNAF  thermochemical  rate  constants. 

For  <fi  = 2.5,  a flame  temperature  of  2250°K  was  calculated, 
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resulting  in  a value  of  p of  6.72  x 10  j sec/m  and  a value 
of  Z of  350  nm.  Results  of  particle  sizing  measurements  are 
shown  in  Fig.  3 for  various  heights  above  the  flat  flame  burner, 
scattering  angle  0,  and  fuel  air  equivalence  ratio  <j>. 

DISCUSSION 

The  results  shown  in  Fig.  3 are  encouraging  because 
they  indicate  that  Diffusion  Broadening  Spectroscopy  can  be 
applied  in  a straightforward  manner  to  determine  particle 
sizes  in  flames.  Measured  particle  diameter  is  observed  to 
increase  from  40  nm  to  260  nm  as  the  height  h above  the  flame 
front  increases  from  0.5  cm  to  3.0  cm.  This  trend  is  con- 
sistent with  the  fact  that  as  h increases,  particle  residence 
time  increases  and  the  probability  of  particle  agglomoration 
increases.  No  definite  effect  of  equivalence  ratio  on  par- 
ticle size  was  observed.  Initial  particle  formation  was 
believed  to  occur  at  a location  h = 0.1  cm,  as  evidenced  by 
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broadband  visible  radiation  in  that  region.  However,  the 
ratio  of  signal  to  flame  emission  was  less  than  10  in  the 
region  h = 0.1  to  0.5  cm.  It  was  therefore  concluded  that  in 
order  to  determine  particle  size  in  the  size  range  below  40  nm 

in  this  flame,  it  is  necessary  to  either  provide  laser  power 

' / 

in  excess  of  28  milliwatts  or  to  reduce  the  bandwidth  of  the 

optical  system  to  less  than  1.6  nm. 

Laser  power  requirements  for  flame  studies  are  dictated 

by  particle  concentration  and  size,  flame  emission  intensity, 

and  the  detection  optics.  As  particle  size  decreases,  the 

intensity  of  Rayleigh  scattered  light  decreases  as  d^,  while 

2 

broadband  flame  emission  only  decreases  as  d , eventually 
necessitating  the  use  of  higher  laser  power  or  improved 
optical  filtering.  The  ratio  of  scattered  signal  to  flame 
emission  can  be  calculated  from: 

2 

! = IP  ~~2~~1~>)  + cos2  6>/Te  RU,T)  d2  AX)  (19) 

where  is  the  laser  power,  AX  is  the  bandwidth  of  the  optical 

system,  Aq  is  the  local  laser  cross  section  area,  m and  e are 

particle  index  of  refraction  and  emissivity,  respectively,  and 

2 

R is  the  Planck  function  in  watts/cm  ym. 

The  effects  of  other  spectral  broadening  mechanisms  were 
assessed  by  varying  the  scattering  angle  0.  Diffusion  broaden- 
ing is  the  only  mechanism  that  has  a 0 dependence,  and  Eq.  (17) 
predicts  that  the  ratio  of  halfwidths  measured  at  0 equal  to 
60°  and  25°,  respectively,  for  a monodisperse  aerosol  is  5.33. 


The  measured  ratio  HW  (60°)/HW  (25°)  was  found  to  be  4.2  at 
location  y = 1 cm  and  <j>  = 2.5.  The  large  0 dependence  mea- 
sured in  this  study  therefore  indicates  that  the  primary 
broadening  mechanism  is  diffusion  broadening.  Time  of  flight 
broadening  was  calculated  to  be  less  than  0.7  Khz.  Further- 
more, the  deviation  of  the  above  result  from  that  predicted 
by  Eq.  (17)  for  a monodisperse  aerosol  is  used  later  to  esti- 
mate the  particle  size  distribution  function. 

An  electron  microscope  was  also  used  to  size  the  carbon 
particulates,  using  a Nichrome  wire  to  obtain  samples  at  3 cm 
above  the  burner.  A typical  electron  micrograph  is  shown  in 
Fig.  4.  Particulates  were  observed  to  be  in  the  200  nm  range, 
in  agreement  with  the  results  of  Fig.  3. 
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SIZE  DISTRIBUTION  MEASUREMENTS 


It  is  reasonable  to  ask  two  questions  concerning  the  use 

of  Diffusion  Broadening  Spectroscopy  in  polydisperse  aerosols. 

First,  what  does  the  mean  diameter  determined  from  Eq.  (17) 

represent?  Secondly,  is  it  possible  to  deduce  the  particle 

size  distribution  parameters?  To  assess  the  effects  of 

polydispersity , consider  an  aerosol  whose  size  distribution 

function  n(d)  is  composed  of  N discrete  size  components. 

The  resulting  diffusion  broadened  homodyne  spectrum  is  the 

2 

arithmetic  sum  of  N simple  Lorentzians,  given  by: 
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Vknink  k2(di  + V 

[K2(Dj  +Dk))2  + f2 


(20) 


where  A depends  on  the  optical  geometry,  and 
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where  a . and  a,  are  the  Mie  scattering  cross  sections  corres- 
] * 

ponding  to  jth  and  kth  particle  diameters,  respectively.  In 

addition  to  N Lorentzians  (one  for  each  size  component) , 

2 

Eq.  (20)  includes  N - N cross  beat  terms  which  result  from 
beating  between  optical  spectra  that  correspond  to  the  jth 
and  kth  size  components. 

To  assess  the  effects  of  polydispersity , a computer 
routine  was  developed  to  evaluate  Eq.  (20) . A log-normal 


particle  size  distribution  was  used,  which  was  character- 


ized by  a geometric  mean  diameter  dg  and  a geometric  standard 
deviation  s.  Mie  scattering  cross  sections  of  carbon  particles 
are  given  by  Fenn  (1966).  Senftleben  and  Benedict 

(1919)  have  measured  the  index  of  refraction  of  carbon  par- 
ticles to  be  m = 1.59  - 0.66i.  Computer  calculations  were 
performed  with  gas  conditions  and  scattering  angles  equal  to 
value  corresponding  to  the  flame  study. 

Results  are  shown  in  Figs.  5,  6,  and  7.  The  measured 
halfwidth  in  Fig.  5 is  strongly  affected  by  the  degree  of 
aerosol  polydispersity.  The  weighted  mean  diameter  dm 
obtained  from  Eq.  (17)  is  therefore  not  a good  estimate  of 

'i. 

geometric  mean  diameter.  However,  Fig.  6 shows  that  the  value 

of  d calculated  from  Eq.  (17)  is  a fairly  good  estimator  of 
m 

the  surface  mean  diameter.  This  result  means  that  a reasonable 

estimate  of  surface  mean  diameter  can  be  determined  in  poly-  * 

I * 

disperse  carbon  aerosols  simply  by  using  Eq.  (17).  A similar  *' 

result  was  obtained  by  Hinds  and  Reist  for  nonabsorptive  
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particles  for  which  m = 1.52.  In  both  cases  the  correlation 
obtained  is  an  empirical  result  and  may  not  be  valid  for 
particles  of  a different  index  of  refraction. 


A method  to  measure  the  size  distribution  parameters  is 
now  proposed.  For  the  case  of  carbon  particles,  our  results 
show  that  measurements  of  homodyne  spectral  halfwidth  at  two 
angles  is  sufficient  to  determine  size  distribution  parameters 
dg  and  s.  Optimum  sensitivity  is  obtained  by  properly  choosing 
the  two  scattering  angles.  In  effect  we  are  combining  the  DBS 
technique  with  well-known  Mie  scattering  principles. 

Figure  7 indicates  the  range  and  limitations  of  size 
distribution  measurements.  The  measured  halfwidth  ratio  is 
seen  to  be  adequately  sensitive  to  variations  in  the  value  of 
s in  the  range  dg  > .025  ym  and  s > 1.5. 

The  poor  sensitivity  of  the  .01  ym  curve  in  Fig.  7 is 

due  to  the  fact  that  the  scattering  cross  section  is  no  longer 

angular  dependent  in  the  Rayleigh  scattering  regime  (d  < X/10) . 

In  the  Rayleigh  regime,  all  of  the  Lorentzians  in  Eq.  (20)  are 

broadened  equally  and  amplified  equally  by  a factor  of  sin  0/2. 

Our  computer  model  verified  that  the  resulting  spectrum  is 

2 

broadened  exactly  by  sin  0/2  and  that  the  degree  of  poly- 
dispersity  cannot  be  deduced  for  d < .01  pm  for  our  conditions. 
Adequate  sensitivity  and  the  single  valued  nature  of  the  curves 
in  Fig.  7 are  due  in  part  to  the  relatively  well  behaved  nature 
of  the  Mie  cross  section  for  carbon;  this  may  not  be  true  for 
nonabsorbing  particles. 
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Values  of  dg  and  s can  be  deduced  directly  fro™  Figs. 

6 and  7.  Measured  values  of  halfwidth  and  halfwidth  ratio 
correspond  to  horizontal  lines  on  each  graph.  Values  of  dg  and 
s are  then  deduced  visually  by  iteration.  For  example,  data 
obtained  at  h = 1 cm  in  the  acetylene/oxygen  flame  used  in 
this  study  corresponds  to  HW  (60°)  = 8.1  KIJz  and  HW  (60°)/ 

HW  (25°)  = 4.2.  The  deduced  value  of  dg  is  .105  + .03  ym 
and  s = 1.5  + .2.  Relatively  large  uncertainty  results  in 
this  case  from  the  double  valued  nature  of  the  curves  in 
Fig.  7 over  a small  range  of  conditions. 

The  research  reported  in  this  paper  was  performed  at  the 
Air  Force  Rocket  Propulsion  Laboratory,  Edwards,  California; 
funds  were  made  available  by  the  Air  Force  Office  of  Scientific 
Research  and  administered  by  Auburn  University  under  the  1977 
USAF-ASEE  Summer  Faculty  Research  Program.  O 
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FIGURE  CAPTIONS 

. Schematic  of  Experimental  Arrangement 

. Measured  Power  Spectrum  in  an  Acetylene-Oxygen  Flame 
(8  = 25°,  H = 3.0  cm) 

. Particle  Diameters  Measured  in  an  Acetylene-Oxygen  Flame 

. Electron  Micrograph  of  Flame  Particulates  Sampled  at 
H = 3.0  cm  (Indicator  Width  = 0.2  pm) 

. Sensitivity  of  Measured  Spectral  Halfwidth  to 
Geometric  Mean  Diameter  for  a Polydisperse  Carbon 
Aerosol 

. Correlation  of  Measured  Diameter  [Eq.  (17) ] 
with  Surface  Mean  Diameter 

. Sensitivity  of  Spectral  Halfwidth  Ratio  to  Geometric 
Standard  Deviation  in  a Polydisperse  Carbon  Aerosol 


